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^ . 

' We find tiiat tlie conditional statistics of temperature difference at fixed values of the locally av- 

' eraged temperature dissipation rate in turbulent convection become Gaussian in the regime where 

(N . the mixing dynamics is expected to be driven by buoyancy. Hence, intermittency of the temperature 

(-H ' fluctuations in this buoyancy-driven regime can be solely attributed to the variation of the locally 

, averaged temperature dissipation rate. We further obtain the functional behavior of these condi- 

• tional temperature structure functions. This functional form demonstrates explicitly the failure of 

\^ ' dimensional agruments and enhances the understanding of the temperature structure functions. 

^ ' In turbulent fluid flows, physical quantities such as velocity, temperature and pressure exhibit seemingly irregular 

a fluctuations both in time and in space. A key issue in turbulence research is to make sense of these fluctuations. 
I The central result of the seminal work of Kolmogorov in 1941 (K41) ||l| is that the fluctuating velocity field in high 
Reynolds number Navier-Stokes turbulence is self-similar at scales within the inertial range, the range of length scales 
that are smaller than those of energy input and larger than those affected directly by molecular dissipation. K41 
^ predicted that the velocity structure functions {[u{x + r) — u{x)]p) scale as r^p with scaling exponents ^j, equal to 
"j^ p/3 when r is within the inertial range. Experimental and numerical results, however, indicate that is a nonlinear 
function of p and that turbulent velocity fluctuations are scale-dependent in that the shape of the probability density 
I , function (PDF) of the velocity difference u{x + r) — u{x) changes with the scale r even when r is within the inertial 
■ O ■ range. This deviation from the K41 results is associated with intermittency or the uneven distribution of turbulent 
activity of the velocity field in time and in space. 

Extensive efforts have been devoted to the understanding of the problem of intermittency or anomalous scaling. 
In his Refined Similarity Hypothesis (RSH) 1^,^, Kolmogorov attributed this intermittent nature of the velocity 
fluctuations to the spatial variations of the energy dissipation rate. Various models have been put forth for the 
' ^ statistics of the locally averaged energy dissipation rate. The most recent model of She and Leveque Q proposed a 
' hierarchical structure for the moments, which leads to predictions that are in good agreement with experiments. This 
t — ' moment hierarchy was later shown to be naturally satisfied by log-Poisson statistics . 

High Rayleigh number convection has been a well-studied model system for investigating turbulence. Fluid motion 
is driven by an applied temperature difference across the top and the bottom plates of a closed experimental cell filled 
with fluid. The temperature field in convection is thus an so-called active scalar. The flow state is characterized by 
^ the geometry of the cell and two dimensionlcss parameters: the Rayleigh number Ra = aglS.L^ j (vk) and the Prandtl 
"Y^ ■ number Pr — v/n, where L is the height of the cell, A is the applied temperature difference, g the acceleration due 
I to gravity, and a, v, and k are respectively the volume expansion coefficient, the kinematic viscosity and the thermal 
^ ■ diffusivity of the fluid. When Ra is large enough, the convection becomes turbulent. 

In turbulent convection, the temperature fluctuations are also intermittent 0. As for velocity fluctuations in 
high Reynolds number Navier-Stokes turbulence, it is of interest to understand the intermittency of temperature 
Q fluctuations in high Rayleigh number convection. Turbulent convection poses additional interesting questions of 
Q its own. There is the issue of whether and how the characteristics of turbulence are affected by the presence of 
^ ^ buoyancy. One expects the mixing dynamics to be driven by buoyancy at scales larger than the Bolgiano scale, 
•'H . = e^/''/[T^'^/^(a(7)'^/^] where e and x are respectively the average energy and temperature (variance) dissipation 

I rates. On the other hand, for length scales smaller than Ib, the mixing dynamics is expected to be driven by the 
^ . inertial force of the fluid motion and the temperature field is effectively passive. Recently, one of us (Ching) |^ 
" " ' has analyzed the intermittency of temperature flcld in turbulent convection. The normalized temperature structure 
functions have indeed been found to have different scaling exponents in the buoyancy-driven and in the inertia-driven 
regimes. 

In our present project, we attempt to understand the intermittency problem of temperature by separating it into 
two parts: the understanding of the conditional statistics of temperature fluctuations at fixed values of the locally 
averaged temperature dissipation rate and the understanding of the statistics of the local temperature dissipation. 
In this paper, we report our study of the first part. The second part of our study is reported elsewhere ||l^. This 
separation allows us to especially address whether RSH type ideas would be fruitful. We shall see that the intermittent 



o 



1 



nature of the temperature fluctuations in the buoyancy-driven regime can indeed be attributed to the variations of 
the locally averaged temperature dissipation rate. Moreover, a change in the statistical features of the temperature 
fluctuations is again observed when the Bolgiano scale Ib is crossed. This change manifests itself as a change in the 
behavior of the conditional PDFs of the temperature difference at fixed value of the locally averaged temperature 
dissipation rate. 

We use temperature data obtained in the well-documented Chicago experiment of low-temperature helium gas 
| pl] , p^ for our analyses. The experimental cell heated from below is cylindrical with a diameter of 20 cm and a height 
of 40 cm. A mean circulating flow is present for Ra > 10*. The temperature at the center of the cell, T(t), was 
measured as a function of time t. We evaluate the temperature difference between two times: Tr{t) = T{t + T) —T{t). 
The intermittency of the temperature ffuctuations is manifested as a change in the shape of the PDF of Tr as r varies. 
In our earlier study of this r-dependence Q], the dissipative and the circulation time scales, Td and Tc, were identified. 
A time scale corresponding io Ib is naturally defined hy tb — TcIb/ L. It was shown that Ib can be written as 



Nu^L 
(RaPr)3 



(1) 



where the Nusselt number (Nu) is the heat flux normalized by that when there was only conduction. Thus, tb can 
be easily evaluated using the measured values of Nu, Ra, and Pr. 

The locally averaged temperature dissipation rate Xr is the spatial average of k|VT|^ over a ball of radius r. We 
estimate it by Xrj which is defined as 
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and can be calculated using the one-point temperature measurements. Here, (u^) is the mean square velocity fluctu- 
ations at the center of the cell. 

We start by investigating the conditional PDF of T,-, at fixed values of Xt- We consider those Tr{t) whose cor- 
responding InxT-(i) assumes a certain value within a small range 6, and calculate the conditional PDFs P{Yr\xT) 
where 



Tr 
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As the conditional mean {Tt\xt) is approximately zero, P{Yr\xT) is standardized with zero mean and unit standard 
deviation. For a given r, P{Yt\xt) is found to be independent of Xt for a range of Xt that contains most of the 
data. The conditional PDFs for different values of t are plotted in Fig. 1. We measure the value of Xr in units of 
X = K{{dT / dt)"^) / {uD . In the limit r — > 0, Xr ~ T^, therefore, the conditional PDF is bimodal for small r, as seen in 
the figure. As r increases, P{Yr\xT) changes from bimodal to a function with one maximum and varies with r. But 
for larger r, it becomes a standardized Gaussian distribution and is thus independent of r. Such a change in behavior 
occurs at T « r^. 

Hence, a change in the statistical features of the temperature fiuctuations is again observed as the Bolgiano scale is 
crossed, demonstrating that buoyancy does have an effect on the characteristics of turbulence in convection. Moreover, 
the physical nature of the presently observed change is clear. We have the interesting result that the temperature 
fluctuations at fixed values of Xt become self-similar and thus non-intermittent in the regime where the mixing 
dynamics is expected to be driven by buoyancy. In other words, intermittency of the temperature fluctuations in this 
buoyancy-driven regime can be solely attributed to the variations of Xr- 

In the remaining of this paper, we shall obtain the functional dependence of the conditional temperature structure 
functions (|7t|^|Xt) on p, r, and Xt- 

It is illuminating to flrst work out what functional form is predicted by simple phenomenology and dimensional 
agruments. One expects Tr, the temperature difference across a scale r, depends on r, Xr, and Ur, the velocity 
difference across the same scale r. In the inertia-driven regime, Ur is related to the locally averaged energy dissipation 
rate by Ur (re^)^^'^ while in the buoyancy-driven regime, Ur is generated by buoyancy: u^/r ^ ctgTr- Hence, we 
have 
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Equation implies that 
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if Tr, Xt, and Cr have the same scahng behavior in r as the corresponding quantities with subscript r in r with 

If the variations of Xt and are both ignored, (j^) imphes that the temperature frequency power spectrum has 
a scahng uj^'s for frequency oj < lob and cj^^/'^ for uj > ujb, where ljb — 2tt/tb- The former scahng behavior was 
reported for the temperature frequency power spectra measured in water and hchum p[ while the latter one was 
reported for that measured in low Pr mercury |]l5| . 

Now we proceed with the analyses. From the result that P{Yt\xt) is independent of Xt, we get 

mnXr) ^ F,{T)anT,Xr) (6) 

where 

a{T,Xr)^VW^) (7) 

By definition, F2{t) = 1. For r > tb, P(Yt-\xt) becomes a standardized Gaussian, thus 

f2P n+1 

Fp{r > tb) = J -r{P^) (8) 
is independent of t. For Td < t < tb, wc find that Fp{T) can be fitted by a power law (see Fig. 2), that is 

Fp{T) « CpT"- Td<T<TB (9) 

This r dependence of Fp echoes that of P{Yt\xt) for r < tb- Using (^), such dependence can be attributed to the 
additional variation of the local energy dissipation rate er even when the local temperature dissipation rate Xt is held 
fixed. The scaling exponents ap are plotted in Fig. 3. Since ao = 0:2 = by definition, ap has to be a nonlinear 
function of p, as is found. 

Next, we analyze the functional dependence of a. We fix r and study its dependence on Xt- When r is not too 
large, a^r^Xr) indeed scales with Xt for a range of Xt that contains most of the data. The scaling exponent &(t), 
however, varies with r. When r is large, the data scatter. Thus, we have 

a{T,Xr)^G{T)xl^^^ (10) 

From the relation Xt ~ in the limit of r ^ 0, one gets 6(t) ^ 1/2 as t — > 0. Indeed, as shown in Fig. 4, 6(t) 
is about 1/2 for t < r^. It then crosses over to an approximately linear function of Inr, and has a value of 2/5 at 
T « Tb- This is, therefore, in contrary to the behavior of a(r, Xr) ^ t^^'^xJ'^ and cr(T, Xr) ~ t'^^^Xjt^ respectively 
in the inertia-driven (r^ < r < tb) and buoyancy-driven regimes {tb < t < t^) that simple phenomenology and 
dimensional agruments would predict [see (|^)]. In Fig. 5, we plot cr{T,XT){XT I x)~^^'^^ for various values of Xt- The 
linear fit of 6(t) in Inr is used for t > Td- The data for different values of Xt collapse to one single curve, thus 
confirming (p^). We take the average of the data to get an estimate of G{t)x'^'^'^\ which is shown in the inset. It can 
be fitted by a power law for t > tb with an exponent about 0.27. 

The temperature structure functions {\Tr\P) are related to the conditional ones at fixed values of Xt as follows: 

poo 

mn = / {\TrnXr)Pr{Xr)dXr (11) 

where Pr{Xr) is the PDF of Xr- Using (|) and (|l^, we thus get 

(|T.r)=F,(r)GP(r)(xf(^)) (12) 

Equation ( p^ ) implies that the change in the scaling exponent of the normalized structure functions {\Tt\p) / {T^)p/^ 
observed, when tb is crossed is the combined effect of the change in the r dependence of Fpij) and (x?'''^''). The 
comparison of (|l2|) with data will be presented elsewhere. 

In summary, we have studied systematically the conditional statistics of the temperature fluctuations at fixed values 
of local temperature dissipation Xt in turbulent convection. We have found that such conditional statistics become 
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self-similar in the buoyancy-driven regime, demonstrating that the intermittency of the temperature field in this 
regime can be attributed solely to the variations of Xt- We have worked out the functional behavior of the conditional 
structure functions (|Tt-|^'|xt)- There is indeed scaling behavior in Xt but the scaling exponent 6(r) depends on 
T, in contrary to what simple phenomenology and dimensional agruments might predict. We emphasize that this 
r dependence demonstrates explicitly the failure of dimensional arguments. Together with the knowledge of the 
statistical properties of Xt, this functional behavior would enable us to better understand the temperature structure 
functions. 
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FIGURE CAPTIONS 



FIG 1. The conditional PDFs P{Yr\xT) versus Yr for Ra = 6.0 x 10^"'^ and Xr/x ~ 0-18 for various values of 
T. r = 8 (dotted line), r = 16 (dashed line), t — 32 (dot-dashed line), r = 64 (circles), t — 128 (squares), 
and T = 256 (triangles). It can be seen that PiYrlxr) becomes a standard Gaussian distribution (solid line) for 
T > tb ~ 70. All times are in units of the sampling time = 1/409.6 s. The conditional PDFs are found to be 
independent of Xt- 

FIG. 2. The logarithm of the normalized conditional temperature structure functions Fp{T) = {\TrY'\xT) / {T^\xtY^'^ 
versus Inr for Ra — 7.3 x 10^" and Xr/x — 0.43 for various values of p. The three time scales r^, tb and Tc are 
approximately 6, 60 and 1750 respectively, and are indicated by the dashed lines. All times are in units of the sampling 
time = 1/320 s. p — 0.5 (circles), p — 1.5 (diamonds), p — 1.75 (triangles), p = 2.25 (crosses), p — 2.5 (squares), and 
p = 2.75 (pluses). For Td < t < tb, Fp{T) can be fitted by a power-law CpT°''' (solid lines) and for r > tb, it becomes 
y/2P /TrT{{p + l)/2) (dot-dashed lines) and is thus independent of r. 

FIG. 3. The scaling exponent ap versus p for Ra = 4.0 x 10^ (circles), Ra — 7.3 x lO"'^" (squares), and Ra — 
6.0 X 10^^ (diamonds). 

FIG. 4. The scaling exponent 6(r) versus Inr for Ra = 4.0 x 10^. The time scales and tb are approximately 8 
and 50 respectively, and are indicated by the dashed lines. All times are in units of the sampling time = 1/160.8 s. 
It can be seen that ^(t) is close to 1/2 for t < Td and can be fitted by a linear function in Inr (solid line) for r > Td- 
Moreover, 6(t) « 2/5 at r w tb- 

FIG. 5. ln(T(T, Xt)(Xt/x) '"'^^ versus Inr for Ra = 7.3 x lO^*' for Xt/x 0.13 (circles), Xt/x — 0.35 (squares), and 
Xt/x ~ 0.96 (triangles). The three sets of data collapse into a single function of r (= G{t)x^^'^'') confirming (|l^). 
The times are in units of the sampling time = 1 /320 s while a is in units of the standard deviation of the temperature 
fiuctuations. Shown in the inset is the average of the three sets of data (solid line), which can be fitted by a power 
law (dot-dashed line) for t > tb (indicated by dashed line). 
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